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ABSTRACT 

The Im'tlaT development of a computer program vjhich implements a 
boundary-detection algorithm to detect edges in digital images is de- 
scribeci. Output products include symbol maps for the original data, for 
the Roberts' gradient of the data, and for the boundaries detected in 
the data. Histogram of the original data values, the data values of 
the bouncjary points, and the data values not at boundary points are an 
output. An initial evaluation of the boundary-detection algorithm was 
conducted to locate boundaries of lakes from LANDSAT-1 imagery. The 
accuracy of the boundary-detection algorithm was determined by comparing 
the area v/ithin boundaries of lakes located using digitized LANDSAT 
imagery with the area of the same lakes planimetered from imagery col- 
lected from an aircraft platform. Relative CPU time costs per data 
point are presented for using the computer program on an IBM 370/145 
with an OS/VS operating system. 
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BQundary~DeteGtion Algorithm for 
Locating Edges in Digital Imagery 

M.J. Russell, D.G. Moore, G.D. Nelson, and V.I. Myers* 

INTRODUCTION 

The sensors aboard the SKYLAB workshop and other space-altitude 
platforms acquire synoptic data for many resource applications. An 
example use is for irrigation scheduling. Since irrigation scheduling 
is based upon field variations, a method to detect the spatial limits 
or edges pf fields is required to conduct a fi el d-by- field analysis 
using digital data. After definition of the field edges, data means 
and variances within fields can be determined. 

A boundary-detection aTgorithm was developed to locate field bound- 
aries. The procedures include: 1) data enhancement by computation of 

Roberts ' gradient values , 2) location of data regions near boundaries by 
applying a threshold to Roberts ' gradient data, 3) mapping of potential 
boundaries in Roberts' gradient data by location of local maximums, and 
4) elimination of local maximums which are associated vn'th noise or minor 
edges by applying a gradient-difference threshold. A preliminary single 
feature boundary-detection algorithm is described and illustrated. 

DESGRIPTION OF THE BOUNDARY 
DETECTION PROCEDURE 

Statement of the Problem 

The.radiance differences among individual ground targets are 
measured and recorded by the SKYLAB S-192 mul tispectral scanner. Many 
of the ground targets are relatively homogeneous clusters of radiance 
values in the data, e.g. a field of corn, a Take, or a fallow field. 
Boundaries between homogeneous clusters are located where abrupt changes 
in magnitudes of the radiance values occur. Therefore, a mathematical 
operator which measures the amount of change in data values would be 
helpful to locate boundaries that separate data classes. 

* Staff Specialist, Remote Sensing Institute; Soil Scientist, Remote 
Sensing Institute; Associate Professor of Electrical Engineering; 
and Director, Remote Sensing Institute; respectively. 
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One Dimensional Model 

An ideal boundary between two data clusters occurs at a step change 
in the data. Ideal data containing no boundaries have small constant 
gradients or have no gradients. Distortion in the form of a point spread 
function and some fandom noise added to this ideal model would form a 
more realistic model of actual data. 

An ideal distortionless and noiseless step boundary along with its 
first and second derivatives is illustrated using solid lines in Figs, la, 
lb, and Ic, respectively {page 11). There is a constant gradient within 
both classes of data and a step change at the boundary. The first deriv- 
ative at the boundary between the two classes is an impulse function and 
elsewhere it is constant. The second derivative is a doublet which changes 
at the boundary from an infinite positive value to zero and then to an 
infinite negative value. The function is zero elsewhere. 

A realistic form for a boundary between two data classes is shown 
with dots in Fig. 1. The infinite slope of the ideal step change becomes 
a gradual smooth change from one class to the other at the boundary. This 
blurring of the ideal data which produces a gradual change at the boundary 
is analogous to low-pass filtering or averaging of the ideal data. The 
first derivative of the realistic data for a boundary between two data 
classes is a maximum value at the boundary between the two classes {dots 
in Fig. lb). The second derivative goes from a large positive value on 
the left side of the boundary to zero at the boundary to a large negative 
value on the right side of the boundary. In summary, first derivatives 
are infinite at ideal step boundaries and are large at boundaries in 
real data. The first derivatives, if no boundaries exist in noiseless 
data, are small constant values or zero. Therefore, a simple decision 
rule for location boundaries in noiseless data is to locate all points in 
the first derivative data that are larger than the surrounding points. 

When realistic data are used for analysis, small changes in data 
values {noise) are normally present within the relatively homogeneous 
data classes. These fluctuations may be caused by actual changes in re- 
flectance or emittance or may be caused by electronic noise introduced 
during data acquisition and/or manipulation. The absolute value of the first 
derivative will have non-zero values at actual data- value fluctuations. 



3 


Boundaries In noiseless first derivative data are located at local 
maximuni first derivative values. A local maximum occurs in the first 
derivative data if a data point has a higher derivative value than adja- 
cent points. The first derivative can reach a maximum not only at major 
boundaries but at locations of noise or minor boundaries. Noise locations 
or minor boundaries would be misclassified as edges using the definition 
that all local-maximum first derivatives are boundaries. The absolute 
value of the first derivative at a local -maximum-point is less for a minor 
boundary or noise than for a major boundary. Therefore, only first deriv- 
atives above a certain level or threshold can be used for initial boundary 
discrimination. Figure 2 illustrates two sets of one dimensional data: 
a boundary in noiseless data and a boundary in noisy data. Maximums one, 
two, and four in Fig. 2b are representative locations of minor boundaries 
located within each of the tv/o relatively homogeneous data classes. The 
smaller maximums within the first derivative data that might be misclassified 
as major boundaries can be discarded. First derivative values for maxi- 
mum one, two, and four are eliminated as major boundaries because they are 
less than the selected gradient threshold indicated in Fig. 2b. 


Two Dimensional Model 


With two-dimensional data, as in one-dimensional data, an abrupt change 
in the magnitude of the data occurs across a boundary. The first partial 
derivative in the direction perpendicular to the boundary measures this 
change, A spatial gradient can be used to measure magnitude changes in 
data. This value is the maximum rate of data change, regardless of direc- 
tion, for a specific data location in the matrix of data. If the point is 
located at a boundary, the direction of maximum change will be perpendicular 
to the boundary. 

For continuous functions F(X,Y) the magnitude of the spatial gradient 
|vF{X,Y)j is define?! as: 




3F 2 
ax 


ay 


The partial derivatives in this expression can be approximated by the 
first differences between data values in two orthogonal directions. A 
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discrete approximatibn of 
defined as: 



this magnitude of a spatial gradierit is 




This discrete form is approximated by using the Roberts' gradient which 
is defined as : 




^i+l,j+l‘ '*'ij+l 


F. 


i+lj 


Abrupt changes iri the original data values which represent boundaries 
become relative maximums in the Roberts' gradient data. Changes due to 
minor variations of data values or noise within the relatively homogeneous 
data classes can also become maximums. However, the magnitude of these 
maximums is usually much smaller than those located near boundaries. A 
gradient threshold eliminates the smaller gradierits that might become max- 
mums. 


The use of a gradient threshold reduces computer time necessary for 
locating boundaries in the gradient values below this threshold. To 
locate a potential vertical boundary, the Roberts' gradient values, which 
are greater than the gradient threshold, are processed to detect if they 


are local maximums in the horizontal Tines. A horizontal boundary is 
located by evaluating the Roberts' gradient values in the vertical lines 
in a similar manner. However, all maximums aren't boundaries. They are 
boundaries only if the local maximums are found while scanning the Roberts' 
gradient values in the gradient vector direction (normal to real boundaries) 
Without knowledge of the gradient vector direction, scanning the Roberts' 
gradient data for local maximums must precede in two orthogonal directions. 
Thus, local maximums are possible both along and perpendicular to the 
gradient vector direction. The perpendicular case generates local maximums 
parallel to the boundary, i.e. within a homogeneous data class. 

A second derivative operation is used to eliminate the local maximums 
which are parallel to the boundary. The magnitude of the second derivative 
will vary as a function of the direction traversed while crossing a 
boundary. In one dimension, the second derivative of an ideal step 
boundary is a doublet. In two dimensions differencing the gradient values 
(second derivative approximation) along the gradient vec.tor direction 
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produces Targe magnitudes on eit^ the local gradient maximum. 

Differeneing the gradient values along directions other than the gradient 
vector direction will result in lower magnitudes on both sides of the 
.Tocal gradient^^^m Therefore » a gradi ent-di fference threshold elim- 

values perpendicular to the 
gradient vector direction (parallel to the boundary). 

In summarys the boundary-de^t^ does the following: 

T) calculates the Roberts 'gradient enhancing the boundaries tn the 
original data, 2) thresholds gradi values to reduce the gradient data 
to only those locations near boundaries, 3} locates the relative maxi mums 
within the Robe gradient values greater than the gradient threshold, 
and 4) eliminates local occurring parallel to actual edges by 
application of a gradlent-differenee threshold. 

USAGE OF THE BOUNDARY-DETECTION ALGORITHM 

The algorithm has been written as a Fortran computer program for 
an IBM 370/145 computer with an OS/VS operating system. The program is 
listed in the Appendix. This section describes data formats, sequences 
of steps, and output products. The effect of the gradient and gradient- 
difference thresholds is illustrated for several choices of each threshold. 

The program accepts positive integer data ranging in value from zero 
to 255. The data can be in either one-byte words or four-byte words. 

The data are read line-by-line from a tape and stored on a disc. Any 
portion of the original matrix of data can be used by specifying the 
starting coordinates and the size of the data block. 

The program supplies a user with 1) a symbol map with the original 
data quantized to 16 levels, 2) a symbol map of the Roberts' gradient 
data, 3) a symbol map of the boundaries, and 4} histograms. The first 
output is a map of the original data values. The 256 possible data values 
are quantized into 16 equal levels, each represented by a different 
symbol. An illustration of this map is provided in Fig. 3. The original 
data was generated by digitizing LANDSAT-1 band 7 (reflected infrared 
-0. 8^1.1 ym) "9 1/2-inch" transparencies via the SADE system at the 
Remote Sensing Institute. The boundary sought was the shoreline of the 
lake. 
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A second output is a map of Roberts' gradient data, The first 
symbol in the symbol list of the program represents the Roberts' gradient 
value which is equal to the gradient threshold chosen. The symbol used 
for a Roberts' gradient value of the gradient threshold plus one is the 
second symbol. This continues until 48 symbols are used. All gradient 
values greater than the gradient threshold plus 47 are labeled by the 
48th symbol. Examples of Roberts' gradient output maps for two different 
levels of gradient thresholds are illustrated in Fig. 4. 

The maps illustrated in Fig. 5 Were the result of applying two 
different gradient-difference thresholds to the remaining gradient values 
after the gradient threshold of 16 was applied (Fig. 4b). Note that the 
boundary in Fig. Bb, is not continuous. As a result* the gradient 
difference threshold was reduced to four. The resultant boundary is 
presented in Fig. 6. The symbols, presented only at boundary locations, 
represent the original data values. For this illustration, the symbols 
relate to film density. 

The final output is a listing of three histograms. One histogram 
prints the number of data points at each value between zero to 255. The 
histogram of original data values at boundary points and the histogram of 
the data excluding those points located at boundaries are printed as 
separate listings. Examples of these outputs are not contained in this 
document. 

The gradient and gradient-difference thresholds must be chosen by 
the user. The number of data points printed on the Roberts' gradient 
symbol map is dependent on the level chosen for the gradient threshold. 

All Roberts' gradient values of the data which are less than the gra- 
dient threshold will be represented as blank areas as illustrated in; 

Fig. 4. The algorithm only operates on gradient values which are greater 
than the gradient threshold. If the gradient threshold is set too low, 
as illustrated in Fig. 4a, more data points must be checked for bound- 
aries than if the threshold were set higher as in Fig. 4b. When the 
proper gradient threshold is selected, only data points near those bound- 
aries of interest will be represented on the map. The ideal choice of the 



gradient threshold is a threshold which eliminates as many points as'pos- 
, sible from the set of potential boundary pointSs but which does not 
eTimihatee^ part of the actual boundaryv ; At this point of development, 
this decision must, bey proyided ^ by the 

computer. ' . 

In summary, the procedure for using this 'algorithm is : C? : [I 

1. Print a map of the original data. 

■ Print, a' map (either^'^ or 

the total matrix) of Roberts' gradient values-- using a 
gradient threshold no larger than 55i of. the potential 
data-value: range and a gfadisnt-difference threshol d at p 

25% of the gradieht threshold.. 

. : y3; ; Determine the- largeSfc^^ 

be .elimi Rated, ,:whi and 

use this value, for the choice of gradient threshold 
for the enti re scene. 

4. Set the gradient-difference threshold at 25% of the 
gradient threshold chosen in Step 3 and_print a new map 
of Roberts' gradient values . ' 

5. If the boundary is not continuous, decrease the value of 

the gradient-difference threshold. ^ 

: ; 6, If the boundary is continuous and contains additional data 
points, increase the gradient-difference threshold. .. . 

Experience with data containing similar contrast ratios and noise levels . 
will provide the user with a "first guess" for setting thresholds which 
are very close to the required threshold values. The data illustrated in 
this report are distributed between one and 255. The gradient threshold 
most often chosen was 16, and the gradient-difference threshold most often 
chosen was four. The boundary map for these threshold values is presented 
in Fig. 6. 
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EVALUATION OF THE ALGORITHM 

Because SKYLAB S-192 data were not available during the development 
phase of this activity, LANDSAT-1 band-7 imagery v/as chosen to evaluate the 
algorithm. On LANDSAT-1 band-7 imagery, a high tonal contrast exists 
between v;ater and the surrounding landscape. Therefore, the evaluation v/as 
completed using this special class of water boundaries. The digital 
data were acquired by digitization of the 1:1,000,000 scale transparencies 
to obtain spatially correct matrices {in contrast to using the digital 
tape products which are not spatially correct). The LANDSAT-1 imagery 
was digitized at approximately 36 data points per millimeter of film. 

This corresponds to approximately 27,7 m on the land surface. These data 
are similar in resolution to the anticipated data products from the S-192 
scanner aboard SKYLAB. A gradient threshold of 16 and gradient-difference 
threshold of four were required to use the algorithm for these data. The 
sequence of maps which resulted after the various steps of the algorithm 
are presented in Figs. 3 through 6. 

The boundary-detection algorithm was evaluated with comparisons of 
water areas measured on the digitized LANDSAT-1 data and areas measured 
by planimetering enlargement prints from Zeiss {12-"inch" focal length) 
photographic imagery taken from a RB-57 aircraft. The boundary determined 
from the digital LANDSAT-1 data was superimposed on Zeiss imagery for one 
of the lakes studied in Fig. 7. This boundary was within 27.7 m or one 
data sampling Interval of the actual boundary. The area of Lake Goldsmith 
of Brookings County, South Dakota, by planimetering from the Zeiss imagery 
v/as 120 hectares of 297 acres. The white line superimposed onto the 
aircraft image is the boundary located, not from the aircraft image, in 
the LANDSAT-1 data. The area determined using the LANDSAT-1 data was 119 
hectares of 295 acres. 

The data values located at the boundary betv/een the lake and land using 
these digital LANDSAT-1 data were distributed between values of 112 to 
192 (Fig. 8). Therefore, a significant improvement over simple thresh- 
olding to locate the v/ater v/as possible with boundary detection because 
no one threshold of data values would have located the same boundary that 
was obtained using t-he boundary-detection algorithm. These results imply 



that for the data used, simple thresholding of the data v/ouTd not locate 
the boundary found between land and water in the digital LANDSAT-1 data. 

EXPECTED ERRORS IN COMPUTING AREAS OF WATER BODIES 
WHEN BOUNDARIES ARE LOCATED USING THIS PROCEDURE 

An accuracy evaluation was conducted by measuring lake areas on 
LANDSAT data using the algorithm and comparing these areas to the actual 
area as measured from aircraft photography. The sampling frequency using 
these input data from the "9 1/2-inch" transparency at 36 points per milli- 
meter corresponds to the ground interval of 27.7 m. For the illustration 
in Fig. 7, the predicted boundaries from LANDSAT-1 data were located within^ 
27.7 m of the actual boundary. Assuming the 27.7 m value to be the max- 
imum expected error for this specific application, the error of measuring 
various sized water bodies can be computed and presented as in Fig. 9. 

The maximum expected error in percentage of total area v/as computed using 
the 27.7 m value for circular lakes. An expected error of less than 10% 
can be anticipated for lakes larger than approximately 3 hectares or 7.4 
acres. The lake used as an example in this report was 120 hectares or 
297 acres and was measured within 0.6%. The are measured error values 
for various sized lakes which were evaluated using similar procedures as 
outl ined. 

A circular lake was the least perimeter per unit area of any geometric 
shape. Therefore, by increasing the perimeter per unit area for irregular- 
shaped lakes, the error in measuring lake area increases. The maximum 
expected error in measuring lakes having differing shoreline development 
factors is illustrated in Fig. 10. 

COMPUTER TIME FOR USING THE ALGORITHM 

The costs for data processing are dependent upon the system used. 

An IBM 370/145 computer with an OS/VS operating system was evaluated. 

The initial compiling time for any size of data matrix was 28.8 sec. The 
additional computing and printing time averaged approximately 1.18 milli- 
second per data point. As an example, the dollar cost incurred for a 100 
by 100 matrix of data (10,000 data points) would be approximately $3.40 
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on the South Dakota State University IBM 370/145 computer with an OS/VS 
operating system. This cost is based upon the initial input data as a 
digital tape. Figure 11 illustrates the computer time requ'^red for ap- 
plication of the computer program for various sizes of data matrices. 
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Fig. 1 - Ideal and distorted data representations for a step boundary 
between two data classes where the solid line represents an 
ideal set of data and the dots represents a distorted data 
set. The ideal data contains a step boundary; whereas, the 
data are gradually changing at the boundary. 



12 


a. Plat of original data values. 

Boundary 



Fig. 2 - Distorted and expected data of a step boundary where the solid 
line is a distorted but noiseless data set and the dots are a 
distorted and noisy data set. 
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Fig. 3 - A symbol map of the original data for a digitized LANDSAT-1 
image. The symbols represent film density values quantized 
into 16 equal levels of data which are distributed over a 
range of zero to 255. The 16 symbols are . ,1 ,2, 3,4,5, 6, 7, 8, 
9, A,B,C,D,E, and F, respectively. They are in order of 

increasing magnitude, where equals 0-15, ,F 

equals 240-255, etc. 
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Fig. 4 - A symbol map of the Roberts' gradient of the digital LANDSAT-1 data (presented in Fig. 
resulting from tv;o different levels of gradient thresholds. The boundary-detection 
algorithm considers only these remaining points when searching for edges. 
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Fig. 5 - The boundary symbol map of the digital LANDSAT-1 data (presented in Fig. 3) with two choices 
of the gradient-difference threshold. The symbols are the same as those in the data output 
map (Fig. 3). The gradient threshold level was 16. 


Note that in 5a the boundary is noisy, indicating that the gradient-difference threshold was 
too small. In 5b, the boundary lost continuity indicating the threshold was too large. 
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The boundary symbol map with the best choice of both gradient 
and gradient-difference thresholds. The gradient threshold 
was 16 and the gradient-difference threshold was four. An 
example interpretation of the boundary location is denoted 
by the solid line. 





Land-lake boundaries on Zeiss imagery ("12-inch" focal 
length camera in RB-57 aircraft) compared vn’th those located 
using digital LANDSAT-1 data. The v/hite line is the boundary 
located using the boundary-detection algorithm applied to 
digitized LANDSAT-1 imagery taken on the same day as the 
aircraft image. 
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Fig. 8 - The boundary located from a digitized LANDSAT-1 image registered 
with a data symbol map. The symbols represent film density 
values quantized into 16 equal levels of data which is 
distributed over a range of zero to 255. Note that a simple 
thresholding of film density values would not have resulted 
in the similar boundary. 
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Fig. 9 - The maximum expected errors for area measurements of circular boundaries if the circle 

radius is overestimate by no more than 27.7 m. The indicate actual error evaluations 
comparing processed LANDSAT-1 data to underflight Zeiss photographic data. The areas 
measured from the scaled aircraft enlargement prints v/ere assumed to be accurate. The 
actual error for these lakes which had a shoreline development factor greater than 1 did 
fall within tolerable limits of the plotted maximum expected error. 
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Fig. 10 - The maximuin expected error for area measurements using 

boundaries located from digitized LANDSAT-1 images of lakes 
with various shoreline development factors. A lake with a 
shoreline development factor of one is circular and a lake 
with a shoreline development factor of five is quite irregularly 
shaped. 
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C SAMPLE MAIN PIOS^AH FCa ThS DETECTlGN OF BCUNOARIcS 

C rRCM DIGITAL DATA 
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c 
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C PURPOSE 

r 

C 01. HEAD THE INPUT MATRIX, WHICH IS OIGITIZEO DATA, 

C FRC'» i TAPE AND WRITE IT ON A DISC. 

C 02.1rl‘l0 THE BG’JNCAP.IES BETWEEN DIFFERENT CLASSES OF DATA. 

C 03-)MND the histogram DF THE INPUT DATA. THE INP-jT DATA v^^•uS 

C DATA POINTS AT BOUNDARl cS ,A hO THE ECUNDARY PCINTS* 
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c 

C.. USAGE 
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C PACKED FORMAT IMPLIES THAT 1 DATA WORD =ONE a PIT BYTE 
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C ! IHE PARAf'ETER I? TELLS THE PROGRAM WHICH.OF TMESt TWO 

C TAPE FORMATS ARE USED) 

C 02. J THE PARAMETER NR, WHICH 13 LISTED ON THE NEXT SHEEI, 

C MUST BE CHANGED ON THE FOLLOWING DIMENSION CARC- 

C t MTECFR-i‘2 NDATCNR-I , I J2 J , IHSGRINR .T RP ) 

r. O’.lTAPr nxrr mu r- >rs •'.“E NE3CEC. {UM T-."Ar fsCOl 
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C flK A8 SYMBOLS. 

C 07.) IF THE DATA IS PACKED SUBROUTINE KSIUN' IS KECESSAPY 
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C DATA. trXflMWLE. .IRc3=l I‘'PLI£S EVEPY LINE =VE«Y MCINT), 

C I tXAM'^LE. .IRcS=3 IMPLIES EVERY THIRD LINE, EVERY U'UD PCINTJ 

C 0B.)ICHAV= if the symbol sequence is TG R£ changed with AVCfHhR 
C DATA card, a one is PLACED HiRE, IF A CHANGE IS T.UT h ANTED, 

C A IS PLACED In THIS FIELD. 


.. Ot.lir-L IF The data is packed, tP = D IF THE DATA iS UNPACKED. 

C 0/.)“X IS T}.E X COOROtNATc CF THE START CF THE DATA BLOCK USEC 
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RDNDOOl 
RdNOOOS 
R3N0010 
RaNDQlS 
RaNUC2E) 
RaND025 
RaN0027 
RBNOCaO 
RBfl003S 
R3N03A0 
SBNOD^S 
RB\0050 
RBMuaSS 
ReNDOBO 
PBND065 
KBNDC70 
R2NDCV5 
R3NC380 
nBSDGSS 
kBNDr'TO 
36NDC95 
RaNGlOD 
fiBNDlO,? 
KG 1 05 

i*o< -u' fc. i. 

RSNUt 15 
RBNDUO 
R3N’'- I2i 
K3NC,;5 
waNC I 30 
RBNCtRZ 
KGND133 
RBf;D134 
hBNO n5 
RiiND 136 
ROND 137 
M3NG1A0 
rt'-NOISO 
RBND160 
SBf.Dl70 
RBN-018 j 
RBNDIMO 
WBNDEOO 
RBND210 
KBNG220 
RBN0230 
RBAJG 2 AQ 
R3N0260 
RRND^TO 
RtlN02b0 
HRNUEqQ 
REN J IDO 
RBNij 320 
RB'|f)32C 
P.8 ■(0330 
RBKT153 



oooor'onrir>firirinr'r'ir»oor>or>ririnrinor>nor>nriooor^oor»<~>nonoort 
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lO.lLRBG^ IS iHc THR = SMC 1 LD FOR THc GRAOIENT VALUE FOR NOISE l<ITHrN 
DATA CLASSES ^C^ NEAR BOUNDARIES. 

A typical range FOR LRBGR... 6 GT LRBGH LT 


... KETHOD 

01. JREAD THE TAPE OF INPUT PA JRI X IP. 

02 . )hHirE THE INPUT MATRIX IH ON THE DISC. 

OJ.IOPCRATE ON THE FIRST 123 COLUMNS OF DATA STORED CM THE DISC, 
123 COLUMNS k»AS CHOSEN DEGAUSS ONLY 128 POINTS CAN BE PRiNTcO 
PER LINE PRINTER PAGE WIDTH. 

OA. IUNPACK THE DATA IF NECESSARY. 

05. JTAXE THL ROBERTS GRADIENT OF THE DATA..* 

ROBERTS GRADIENT =lABSlNDATtI,J)-NDAT(H-l,J<-n) + 
IABS(NDArn.J*l)-NOAT<I«-l»JJ ) 

06 . )-.RITE A PAGE OF DATA NDAT IN SYMSCL NOTATION. 

UUANTIIED INTO 16 LEVELS. 

07 . )wRITE A page of ROBERTS GRADIENT !N SYMBOL NOrATION . 

SYHllMIS equivalent TO LRBGR 
SVMin= IS EQUIVALENT TO LRCGRi^i -1 

OB. ITHE BOUNDARY IS FOUND FP-OP RCBFRTS GRADIENT INFORKATIQN 

DATA IN GOTH THE X AND Y 0 ! RECTI DNS. 

O'i.IA BOUNDARY CANNOT BE FOUND FRCM ROBERTS GPADl ENT • I NFQRPAT lON- 
IRBURIUJ) IF iHc gradient IS LESS THAN THE THRESHOLD LRBGR. 
lO.JwHEN A MAXIMUM IS FOUND WHILE WORKING IN THE X OR Y DIRECTIONS 
THE MAXIMUM DIFFSaENCE IG IS CALCULATED BETWEEN THE PAXIHUH 
AND ThL GRADIENT POINTS ON EITHER SIDE OF THIS POINT 
IN THE X OR Y DIRECTION RESPECTIVELY. 
ll.ITHE ThRESHOLD NTHS ELIMINATES ALL BOUNDARY POINTS IN WHICH THE 
i-iAXiMU” DIFFERENCE IG AS CALCULATED ABOVE IS LESS THAN NTHS. 
L 2 . 1 THE BOUNDARY IS CLASSiriEH BY FINDING THE AVERi-^E CF THE 
HliUR uaTR PiliNTS <iu Tnt Etu:'>GART FCINT. 

13 . )A SYMBOL PRINTOUT QUANTIZED INTO 16 SYMBOLS CF THE DATA 
POINTS AT A boundary IS PHISTtD ON THE lI NE PRINTER. 
lA.JTHE ABOVE PROCESS IS REPEATED IF NECESSARY ON THE T tCOND 12 R 
COLUMNS OF DATA. 

I 5 .JTHE HISTOGRAM OF THE DATA, THE HiSfDGRAM OF THE DATA MINUS 
THE BOUNDARY POINTS, AND THE HISTOGRAM OF THE DATA IS 
NOW PRINTED OUT, 


INTEGCR-OZ IP.BGRIN0,132 ) , NDA T ( <30 , 1 32 1 
I.NTEGER’a2 NA 1 1 32 1 , NQ ( 1 32 J 
DIMENSION SYMI48) 

CIHENSION DEC! 154» 
r.IHENSION KCCUNH256; 

DI'^ENSION JC0UNTI236J 
DIM-NSID'. IC0UNTt256) 

CQM'TN/COMnrtC/LlNEI 2C'.8 3 /COMB YT /L I NEPK ( 5 12 J 

LATA SY*'/' 11 n '2222 ' t • 3333* . , *5555 • , 

L ' t icL* , ' 17 77' , ‘fleSB’ , •RJJqP • , »AAAA* , » “23? * . ‘CCrT* , • CCCO* , 'FEEE* , 
d'FF f f * , 'Gf.:U', 'HPHH* ,'lin*,'JjJJ»,*KXK%',*LLLL',* •'PM,'" .'NNNN' , 
6'fir J M , 'Pppp ' , *OCCQ •,*RR.>R*,»S5S£*.*TrTT'. HICUU* , »>*VVV* , 'WNUW* , 

6*xxAX*,'YYVy','ZZzz*,’,,,,','‘***',* ' .'* = = = •, 'Vt***,*!!!!*, 

'CCCO *f. Etc, 'MS**/ 

IC-(-= CA'-C reader as INPUT. 

LP= LINS PRINTER AS OL'TP'jT. 

p.'uNCTir tape as ippur. 

IC‘ ^ i . 


RBND335 
RBND 3',0 
RBND 3'»3 
RbND2A5 
R6N3 JAB 
RBN03S0 
RBND355 
RBMD360 
RBND370 
RBn 0380 
RBN0390 
HBNDACO 
RBND610 
BBN0430 
RBNDA40 
R3ND650 
RBN0A6G 
RBNU470 
RBNDA80 
ReND-igG 
RBND500 
RBND51C 
RBND530 
RBND540 
RBND550 
R8N056D 
ROND 5 70 

RarjosBO 
RBND590 
RBND600 
RBN0610 
bND6 1 S 

' . j ,1 M L kJ 

R3NC620 
RB.\Dt.30 
tt 3 K D £3 *♦ Q 
H 3 ND 6 SO 
RBfjD660 
R5ND665 
•■’,BND670 
R3N0680 
,RS'I0690 
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LP=12 

MAGT^O'J 

S-;tL I IC.’’ f 9 )LCL .NK,\2X, IRES , I CHAN ,I P , “X , r' Y , N f HS , L RBGR 

IF ( 1 CHAN.EQ. n PEA'J ( ICa , <,2) SYM 

LBL>'-L3L/<. 

NR =\R- I Kc S 
NRK^NR/ IRES 
CD 1 IH=1,?55 
KCO'JNT ( I h) =0 
1 ICCiJMI IH) =0 
RcWINO 8 
LQ S13 1=1, MY 

bJ3 REARI'-AGT.^O ) I L I \E PK( I 0 » . ID = I , L8LH) 

C RE4D TAPE v,3 I TE TD DISC 
CO 0 3 1 = 1, NR, IRES 
CO ll IG=1,1RES 

11 READI^AGT, AO ) ILINEPK.I ID) , ID=l ,L8LWI 
0 3 hRI T£( 3) (L INEPKl 10) , I 0= I , LBL V, ) 

CO 8 LM^1,NEX 
V.R1 irtLP, IC) 

L = HX*1 130* IRES) 

I F( Lf . GT . 1 . ANO.L .LT. L8L) L = LM 

IF(L^'.EO.NEX.ANn.lP.EO.O)L=LflLW 

IFU v.EO ,M Ex. AND. I? .£0. l)L=L8L 

JC= I U-MX ) / IRES *l 

JCC=JB-1 

J<U= JCU 

IP( JC0.GT.128) JK0=128 

R E « I NO d 
-PI rt(l.P,3CA)LM 
CO ** J iJ = l,NRR 
C TucnAIA 

R'. AOl fi ML riCPKI ID) , I 0 = 1 ,L3LW) 

IM IP.EO. I )GU lO 6 
LC ! -Ki)=I,LBLW 
7 LINEI-KG)=L INEPKI JKO) 

CO TO 12 
S CALL r.ilONP 

12 IFUJ.t.EO. 1 ) GO TO 600 
CO 21 KK=l,JD 

21 NAC- K )=Ni3IKK) 

800 PB=0 

CO S !.K = MX,L ,IRES 
PB=Md .1 

NOAM JJJ.MP)=LINEIK;;) 

5 .\6(>I«)=L INEIKX) 

C P-D-JEKTS GRADIENT 

I FIJJ J EO. 1 ) GO 10 A 

CO 22 l-.KK^UJCD 
JE=KK<Cl 

I YkT s’l-U XK K 1 -N8I JF T 
IXRT =NA{ j£ )-N3UXK) 

A= : Af.:( I tRT ) C! AUil I XKT) 

22 IRftl.A! JJJ,KKK)=A 

1. --irf l!*LDAIA^lOAT 

lO "1 N' = I . JRO 
N’iT':. '.X ) 

Jf.-l ..NT',/ I '■) *2 
IF ('‘i(jT.,.».w .C)J0 = 1 
GECIf )»SYV( JO) 

ICO, •.(<’; .lOMr.TLNT :%JTS5Cl 
t» .t.TS J = :C2./.T C.'JTSI 



H 


(, telllH G^ADl&Nf 

hrtl TEUP. lot 
WHflEttPtlOSJLM 
CD ^5 l*2iN^R 
CQ 2fc J-ltJKG 

ItJ}^LRBGR 
IFUE.ivc.'9ijb»*8 
IF( JD.L1 . UJB-1 
2b DEGI J » = S¥P Ub) 

2 5 «R I T E ( L f , 2 ) ( DEC H Mi) , I N *l , JCO ) 

MiisIP cS-»LM*l2j?l 

jge=jc6-i 

NRE = N'tA-^l 

hfl I rliLP. 10) 

KftlTEClPilO) 

^.KltittP. 306)LM 
CO 29 
|PLUS»l*l 
IM INS* 1-1 
E0 27 J=2,JCE 
JPL US^JC I 

iFQ»IRBGR( 1 1 J) 

IGrO 

II-O 

Jb=l 

IF( IFD.LE.iReaO) GO to 256 
£ CFSGn FO-< X BOU'^GAftllS 

5 1 IPi|»BG<llPl.US.JI 
19 5= M BG3 ) iMiN St J) 

!KM5.;;T.iFOvOa.inP*r.t. |iFO» go fO 5? 

IG =IF0 ^tBU 
Irt^lFO-lH 

IFItP.GT.IO IG>IB 
5 2 CDf.t IMUL- 

C CI4£CK FD« V nOUNCAAIES 
fi3 10> IHJGHI I, JPLUS) 

100> I'RDGRd iJMlNS) 

IF ( I C .GT . I Fd .OR . I SO.Gf ^ If 6 1 CO tO 

llalFD-^tOD 

ID* I FO- ID 

IFMD.Gt.I I ) II*ID 
5t CONflNUE 

I'Ftn.Gt.lC) ]G«M 
IF) ;G.LE."JTHS)G0 to 2S0 

J C* MD& r { I , JM I N s I «NDA T ( 1 1 J ) «NOA T 1 1 PLUS • JH IMS ) > MO A 1 1 1 PLUS . J ) 
X = FL3AM JCI/A.O 
JD* I F IX) XI 

JCCUNT « J D ) - JCCUJN T 1 J O ) ^ I 
KCOJMTUDI =KCDUNTIJ0)*1 

X*X/1.5*6 Cl 

JO= If IX tX) 

IFt JD.LE. U JO*l 
25Q CECIJ )*SVMt JOI 
21 CONTIMUE 

k + n 5t £► , 2) (GEC ( IP) t lM*2i JCEI 
29 CCF.T in JE 
P C-jMl'JUl: 

c fl 'h i'isrns.tftM of the data noAr 

I itiLP. lei 
I I dLi* * 30 31 
W“ J TtJL.'d 201 I 

U*« ! ! {{.?', .€?! 



z • >,12^41) 

'1 f Er .at ; 1(JI5) 

i-:> Ffirti'.it* • i ■ ) 

^ vt f;i;Vv U t • q • , 5 X , 4 1 » S I 5 1 
4.0 F fi-< - ^ r ( 3 ( 2E 0 4 h > ) 

^2 FQ« 'A M 

1 f( • c* ,5 a, Ai ,1 5J 

13£i Fin* AI( * C* , |4.,lAX, 16, 10«, t fi,5X ,16) 

ici Fi,'<-"iri* ■ , IXifcHVEHfl£R,2* ,1 7riTIf'>S SY^SCl USEOl 

302 FQK.'- Ari’ ' , 2 iX, ’QATA' , 3 Xt*i)'AT 4 I'lVUS BGONE AX Y* ,SX, • BOUNyARV* ) 

331 ^wV-.-ATi *0' ,3CA,‘0AfA HlitMSAM'i 

;oi rpi3'£in 31 2CG4A) ) 

336 FaxvArt* •,6CX»*aATA PASE’,12) 

336 F'JJiATi* ',5GX, ’KdecKTS 0*4 Ai) IE NT PAG£*.12» 

136 HIXMATI* • ,sex, "B0»J.NEA.-<Y data P toe ’,12) 
iE3e FQ-XMAT! 3i2eeA6n 
stop 
END 


* SUBROUTINE RSIUNP 

« TWO AP.KAVS AXE INVdLVEp IN THIS SUBROutlNE EACH (vF Ui.tCH StDPES 

-• 2066 eHARAeTERS, 'CQMHKO* STGRES ONE CHARAOtER PER hORD 1 RIGHT 

* JUST IF I EG) t*.‘D •C0MOVT* STGReS UNE PER BYTE. 

* tmS SVICROUTIVc FliRSf ZEROS CGMWRd, THEN Hfl'/ES EACM BYfE ftiQf 

* GE1‘3YT TO THt £CR« E SPOND I NS KCHO IN GDPxRD AND STORE' it KIiGHT 

* JUSTIFltD. TmJS IT UNPACKS THE DATA FrOP ThE flVtE ARRAY AND 

* STORES IT RIGHT JUSTIFIES IN THE U0RC ARiRAY. 

* ENTERiO 3Y 'CALL RSIONP* 

* calling PRQGRA4* KcQ J IKES COPf'EN STATE **EMT 

'= G0PMO : /COM w K 0 / L 1 NE ( 206 a ) /£ G VT /L I N6 RK I 512 ) 

* »EGy„so» and *CCy3YT' ARE sEeUlREb NANES 

6 'ilNt* A.-jD 'LlN^PK' ARE ARtITGarY BUT 

* » 2 P 68 * AND ♦ 612 ' ARE AceUjKttl 

* pROGFAN IS aeusAaiE 

* REMUrt SENSING INSTltUIE 


* SG'JTK DAKCTA StATt UNtMEnNlf .' 

♦ BHOeXINGS SOUTH D Ah OTA 5700e 





pRri'f 

NGGEN 


ENTRY 

R S I UNP 


EXTRN 

GOMkRO 


EXMN 

eOflBYT 


USING 

• ,IS 

asiuNP 

STM 

16, 12. 12 It 3) 


L 

5,=»A<C(jMV,R0) 


L 

7,»A (eUPBYT) 


L 

3,=F*20'*6* 


L 

6 , » F • 6 • 


L 

6 , =F • 3 • 


L 

a,*F*i» 

7 ERO 

fVG 

0 1 6 , 5 ) , «F ‘ 0 • 


AK 

5,6 


PCI 

3, ZERO 


L 

3,iF*2068* 


L 

5,^A<GG.iv,aD> 


AR 

5,6 

MOVE 

IVC 

on .5) ,G( 7) 


ti G T 

3,NEj T 


LH ‘ 

2, 12.281 m 


PV ( 

m 13) ,X'FF* 


DR 

16 

ti E / i 

A« 

S,6 


AR 

7,3 


1 

•. ■; ] 

N* ‘ -t 
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R15 BASE REG. f TER 

SAVE FEXtRAN XEGISTgr 

PUT AP3RESS OF FIRST XGHU IN HS 

PUT ADDRESS CF Fl tST 3Y I £ IN R7 

PUT 206« IN R3 

PUT 6 I*N R6 

PUT 3 IN R6 

PUT I IN Rfl 

• xErg 

WORD 

ARRAY • 

PUT 2068 IN R3 

PUT aG DRESS GF FIRST WORD IN R5 
AD0£=5S LAST BYTE HF THE WORD 
KCVi FREP BYTE ARRAY TC- WQP.D ARRAY 
HAV£ ALU OAIA StEN GUNVcRTED 

♦ IF SO PE TURN 

G0Nt<EL TO 

FOR TP AN ORE 6 ? AN • 

PUT ATriRESS CF NeXT EvTE If; «7 


imODUCKUiufy 



